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Abstract 

In this paper we have demonstrated the double Q 2 -rescaling model (DQ 2 RM) 
of parton distribution functions of nucleon bounded in nucleus. With different 
x-region of I- A deep inelastic scattering process we take different approach: 
in high x-region (0.1 < x < 0.7) we use the distorted QCD vacuum model 
which resulted from topologically multi-connected domain vacuum structure 
of nucleus; in low x-region (10 -4 < x < 10" 3 ) we adopt the Glauber (Mueller) 
multi-scattering formula for gluon coherently rescattering in nucleus. From 
these two approach we justified the rescaling parton distribution functions in 
bound nucleon are in agreement well with those we got from DQ 2 RM, thus 
the validity for this phenomenologically model are demonstrated. 

PACS number(s): 12.39.Ba, 12.90.+b, 24.85.+p 
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I. INTRODUCTION 



Deeply understanding the parton distribution functions (PDF) of the nucleon bounded in 
nucleus are very important both for studying particle and high-energy nuclear physics. Since 
the middle of eighties many important progress has been achieved both from experimental 
and theoretical studies 0. It has revealed experimentally that the PDF of a bound nucleon 
has remarkable difference from that of a free one: the EMC effect in high x-region (0.1 < 
x < 0.7) [0 and the nuclear shadowing effect H] in low x-region (x < 10 2 ) are the most 
important examples. Of the EMC effect, many models have been proposed [0,0, all of them 
could more or less explain the relevant data, but for the nuclear shadowing effect only the 
constituent quark Q 2 -rescaling model and the extended x-rescaling model can make out 
satisfactory explanations after introducing some shadowing parameters As far as we 

know none of these models could describe the PDF of a bound nucleon which cover both 
the available kinetic x and Q 2 range and could explained well all the data from high energy 
I- A and h-A basic processes. Of course building such a model is very difficult, but is also 
very needed. We have tried to take this task and proposed from phenomenology the double 
Q 2 -rescaling model (DQ 2 RS) |J, from which the calculated results agree rather well with 
experimentally relevant data. 

In this paper we want to demonstrate our DQ 2 RM from physical observations. We take 
different approach at different x-region of I- A DIS. In high x- regions (according EMC effect 
which means (0.1 < x < 0.7), for fixed Q 2 , it should correspond relatively low energy region 
so the non-perturbative QCD effects would be important in this x-region. We observe the 
QCD vacuum around the bound nucleon is modified topologically by multi-connected domain 
vacuum structure of the nucleus. We further use soliton bag model of Lee et al |7[] and find 
the radius of a nucleon bounded in nucleus with nucleon number A is larger than that of a 
free one, > rjv 0. According to uncertainty principle, the PDF of a bound nucleon would 
then be soften in comparison with that of a free one and shift the PDF to small x side, which 
are consistent with those from DQ 2 RM, thus we demonstrated our model at high x-region. 
In low x-region which covered by available data (10 -4 < x < 10~ 3 ), it is well known the gluon 
distribution predominated in PDF of a nucleus and the coherent interaction of gluons from 
neighbour nucleons with same impact parameter in a nucleus becomes important. Following 
the method of Levin et al |J , we adopt the Glauber (Mueller) JTOj] multi-scattering formula 
for calculating contributions of gluon coherently rescattering in nucleus, thus we get the 
gluon distribution function of the nucleon bounded in any nucleus (A) and find it agree well 
with those from DQ 2 RM. 
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The organic of this paper is as follows: in sec. II we recapitulate the DQ 2 RM and its 
relatives. A detail demonstration on DQ 2 RM both for high and low x-region will be given 
respectively in sec. Ill and sec. IV. The last sec V is the conclusion and some comments. 

II. RECAPITULATION OF DQ 2 RM 

Since nuclei are formed from weakly bounded nucleons, one could expected the basic 
properties of the nucleon in nucleus, including the behavior of parton distribution, are not 
far from that of the free nucleon. All nuclear effects revealed in I- A DIS should be attributed 
to this nuclear environment which comes from the weakly bounding mechanism. According 
to this observation and informations from relevant phenomenology, we proposed in || the 
following relations of PDF between bounded and free nucleon: 

<lv i (x,0 2 ) = <lv i ( x >Zv<?), (i = u,d,s) 
qt(x,Q 2 ) = qi(x,Q 2 ) = q$ i (x,t£Q 2 ), 

g A (x,Q 2 ) = g N (x,^Q 2 ), (1) 

where rescaling parameter's (a = V, S, g) are depending on and very slowly increased with 
nucleon number A. 

Since it must respect total momentum conservation, 

J xdx S^lQv^Q 2 ) +9 A &Q 2 )} = !■ (2) 

Only two among £, A (a = V,S,g) are independent, so this model is called the double Q 2 - 
rescaling model. 

At present the statistics and the precision of experimental data on relevant processes 
are insufficient to test very precisely theoretical calculation, so in our calculation only the 
leading order approximation on the formulas of hard subprocesses and parton distribution 
functions are considered. Thus in DQ 2 RM the structure function and gluon distribution 
function of bound nucleon respectively is 

F 2 A (x, Q 2 ) = £ e 2 x {q»(x, tfQ 2 ) + q»(x, ^Q 2 ) + qg(x, ^Q 2 )} (3) 

i 

xg A (x,Q 2 ) = xg N (x,^Q 2 ). (4) 

F A (x, Q 2 ) denotes the average nucleon structure function of an ideal nucleus with N = Z = 
\A. 
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The nuclear effect of nucleoli structure functioned and the gluon distribution may be 
represented by following ratios 



R A ' D (x,Q 2 )=F A (x,Q 2 )/F 2 D (x } Q 2 ) } (5) 

R A I D (x, Q 2 ) = g A (x, Q 2 )/g D (x, Q 2 ), (6) 
For the Drell-Yan process in p-A collision the nuclear effect is represented by the ratio 

TAIN = i ^-y x „Q*) i ft^fl 

J dxdxt J dxdxt 



where 



Three parameters £^ (a = ^ ^ d) 1X1 DQ 2 RM are determined as follows: firstly we choose 
suitable and £ A to explain the nuclear effect in I- A DIS process Eq.(|]); secondly we 
describe p-A Drell-Yan process with Eq.(^) and finally determine £ A by using nuclear mo- 
mentum conservation Eq.(f2|) and explain the nuclear effect of 7*-A electroproduction J/ip 
process in Eq.(|B|). 

We adopted the PDF of free nucleon given by Gluck- Reya- Vogt (GRV) JTTJ because this 
formalism allows to be evoluted inversely to Q 2 < 1 GeV 2 , this is quite important point for 
getting the PDF of sea quark and gluon of bound nucleon by Eq.(|l]) at £ A < 1 and £ A < 1. 

The Q 2 -rescaling parameters £ A , £g and £ A on nuclei C 12 , Ca 40 , Fe 56 and Sn 119 obtained 
by above method from the data of references || 12-14] are listed in Table lQ. It is worth to 



note that all of £y and large than unity and increased slowly with A, but all of £ A and £ A 
are smaller than unity and decreased slowly as A increasing. 

Using the parameter listed in Table 1 and relevant Eqs. (l)-(8), together 
with GRV parametrization for PDF of free nucleon, the theoretical predications on 
R A / D {x,Q 2 ),T A / N (x,Q 2 ) and R^ D {x,Q 2 ) for nuclei C 12 , Ca 40 , Fe 56 and Sn 119 by DQ 2 RM 
could be obtained. 



^In ref. || the Q 2 -rescaling constants of Sn 119 which we have fitted are £y n = 1.57, £l| n = 0.45 
and £y n = 0.58, but after recalculating we find it is better to take these values shown in Table | 
instead. 
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The results calculated from DQ 2 RM in ref. || as follows: 1. The curves of R A / D (x.Q 2 ) 
on l-A DIS processes for nuclei C 12 , Ca 40 , Fe 5e and Sn 119 within Q 2 range given in are 
separatedly plotted in Fig.l(a-d) where corresponding experimental data @JT^] are shown. 
2. The curves of T A ^ D (x, Q 2 ) on p-A Drell-Yan process for nuclei C 12 , Ca m and Fe 56 within 
region of 0.025 < x < 0.30, 4 < mg < 9 GeV and Eqm = 40 GeV are separately plotted 
in Fig.2(a-c), where corresponding experimental data are shown. 3. The curves of 
Rg n / C (x, Q 2 ) (the ratio of gluon distribution function of nucleus Sn to that of nucleus C) is 
plotted in Fig. 3 where the experimental data [|14[] of the ratio measured in 7* + A — > J/ip + X 
process are shown. 

For comparison we have also plotted in Fig. 1-3 the corresponding curves which calculated 
from original single Q 2 -rescaling model (as pointed by Li et al in [[HJ. This model with a 
common Q 2 -rescaling parameter £(A) for valance quark, sea quark and gluon, leads to non- 
conservation of nuclear momentum). 

From above illustration, we can see that by using DQ 2 RM all calculated ratios, which 
representing the nuclear effects on PDF of bound nucleon, agree rather well with that from 
relevant experimental data, Since the essential elements of this model are embodied in Table 
1. we shall demonstrate it in following sections. 

III. DEMONSTRATION FOR DQ 2 RM IN HIGH X-REGION 

In Z-A(N) DIS process the Bjorken variable x = §+q2 > where S is the square of center- 
mass energy of 7* and scattered nucleon. For fixed Q 2 , the higher x region correspond to 
lower energy region. Thus in high x-region (0.1 < x < 0.7) the non-perturbative QCD effects 
should be important and the behavior of PDF of the nucleon should depends on space-time 
property of hadron environment. In the following we will use the soliton bag model of T.D. 
Lee [0 to justify the QCD vacuum structure of a bound nucleon is topologically different 
from that of a free one (an isolated nucleon). The consequence of this difference is the radius 
of bound nucleon would be large than that of free one. i.e. ta > r/v- Thus, in comparing 
with free nucleon, the PDF of bound nucleon shift to its small x end. 

comparing QCD with QED and introducing the concept of colour dielectric constant k, 
Lee found that the vacuum polarization of QCD shows a colour anti-screening property. On 
account of this a phenomenological scalar field cr(x) is introduced into Lagrangian L of Lee's 
model to describe the colour dielectric property of vacuum. This a(x) is treated as a classical 
field with an effective potential shown in Fig.(|]). The potential energy V(a) has an absolute 
minimum V(a vac ) = at a = o vac . which represent the QCD physical vacuum out of the 
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nucleoli bag and which correspond to a complete colour antiscreening with k = 0. Inside the 
bag there are QCD perturbative vacuum around the colour quarks, where the field a = 
and correspond to k — 1. 

In the case of nucleus, since many nucleons densely populated, the physical QCD vacuum 
out of any nucleon which bagged by the nucleus is a multi-connected domain and so it has a 
different topologically vacuum structure comparing with that of free one (isolated nucleon) 
which is a single-connected domain. In terms of Lee model which means, for a nucleon with 
A nucleons, the a field in the domain occupied by the remainder (A-l) bags has a value 
a = 0(k = 1) instead of o = cr vac {^ = 0). On the average, the effective QCD physical 
vacuum out of bound nucleon bag correspond to scalar field a a ( the relevant vacuum colour 
dielectric constant ka) which must smaller than a vac , the difference (a vac — a a) increase with 
increasing of A. From Fig.4. one could see, comparing with isolated bag where the vacuum 
pressure on its surface is B and its equilibrium radius is r N , the bound nucleon bags would 
bear smaller vacuum pressure Ba {Ba < B) and enlarge equilibrium bag radius ta{ta > r N ). 
To estimate ABa, a liner relation between a and k assumed as in Lee model, 

= cr vac {l - k a ) (9) 

expand V(cta) closely around a vac and take the first non-zero term, 

V(a A ) = b(a A - o vac f (10) 
where b > according to the definition of physical vacuum. Expressing a a by Eq. @ we get 

AB A = ba 2 vac K 2 A = Gk\ (11) 

A simple and intuitive method estimating the effective constant ka is to consider the 
volume ratio which the nucleus occupied with that of (A-l) surrounding nucleons. From 
geometry picture consideration we get 

A a -1 

Airp A r\ 

where pa is the nucleon number density in the nucleus A. In order to avoid overestimating 
the effect of those bags far away from the considered nucleon a parameter a < 1 is introduced 
and we find | < a < g are reasonable values. 

The relation between ABa and ta is determined by Eqs. (|TTD and (12). On the other 
hand the equilibrium condition of the pressures inside and outside of the bag is determined 



by the formula of the bag model Wfi 
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M = j + yr 3 B (13) 

where c = 3e — a c Z)(o"j£7j)(Af A|) fi — Z is constant in our case, from Eq. (|l3|) we can induce 
to another relation between ABa and ta'- 

3M N r N/ 1 



A5 A = — ^(— --) (14) 
lo7r rt r r. 



Mat is the mass of a free nucleon. Combining Eqs.([TT|), (|12|) and (|14"D we can obtain the 



self-consistent equation for determining the nucleon bag radius, ta, in the nucleus; 

^-i 1 (15) 

since G = fecr^ ac > 0, Eq.(^) shows r A > r N . 

After the EMC effect discovered experimentally, Jaffe and Close et al Q in their Q 2 -scale 
change model proposed that the effective single nucleon structure functions in a nucleus (A) 
are related with that in deuteron by 

F 2 A (x,Q 2 )=F 2 D (x,UQ 2 ) (16) 

where £a is a scale parameter weakly (logarithmically) dependence on Q 2 . On the basis of 
bag model and QCD they further argued that £a is related to ta in nucleus by 

£ A (Q 2 ) = (l!) 2fa (Q 2 /A 2 )/MM 2 /A 2 ). (17) 

By solving Eq.(|T5D we get the radius ta for bound nucleon bag. There are three param- 
eters rjy, G and a in Eq.([T5|) which can be determined through fitting the £a values of five 
nuclei He, Be, Al, Fe and Au, which are obtained by rescaling experimental data; thus we 
get r N = 0.865 fm, G = 0.0362 GeV fm, a = 0.541. In calculation we used fi 2 = 0.63 GeV 2 , 
A = 0.250 GeV. The results calculated for Q 2 = 2.5 GeV 2 are summarized in Table [TT]. 
The experimental values of £a are not so definite. The difference between ^-calculated 
calculated and ^-experimental data for any nucleus A is around 5%. 

Eq. (|i~6|) is a satisfying phenomenologically relation which all the data for different nucleus 
on l-A DIS process are in accordance with it. 

Comparing the parameters = V, S,g) in Table [I] which comes from DQ 2 RM with 

above the parameters £a which describe the effect of QCD vacuum distortion in structure 
function of bound nucleon, we notice that: 
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1. Both the values of in Table | and £ A (in Table ||) are larger than unity and are 
weakly increasing as the nucleon number A increased. 

2. For any nucleus (A fixed), always larger than £ A , the differnece between them are 
essentially unvaried with A. 

Furthermore, after taking the and £g from Table | and using the PDF of leading order of 
GRV parametrization, from Eq.(^) we get F A (x, Q 2 ) for various nucleus. Use Eq.flTSD once 
more we get again the Q 2 -rescaling parameters £,pn2 RM for structure functions of bound 
nucleon from DQ 2 RM. We find for any nucleons, the difference between Cdq 2 _rm with £ A in 
Table are less than 5%. This is just because both the parameter £<j and £ A in DQ 2 RM are 
less than unity and they are decreasing as A increased, as seen in Table fl 
Thus completed our demonstration for the DQ 2 RM in high x-region. 

IV. DEMONSTRATION FOR DQ 2 RM IN LOW X-REGION 

From any phenomenologically parameterization scheme for PDF of a nucleon we know 
that in low x-region, e.g. x < 10~ 3 or even smaller, the distribution function of valence 
quark almost equal zero whereas that of the gluon and sea quark are increased very rapidly 
with the decrease of x; in addition the distribution function of gluon will nearly be 3 to 5 
times as that of sea quark so as first approximation we will neglect their contributions. 

In discussion of high energy 7*-N(A) DIS process at low x-region, it seems more intuitive 
and convenient to use the space-time description approach in target rest frame, Pn(A) — 0, 
rather than as usually use the Feynman-Bjorken infinite momentum Pn(A) oo approach. 
In Pma) = frame, the incoming virtual photon state 7* at a given instant of time, can be 
expanded in terms of its (bare) Fock states [ jlTf : 

|7* > P hys= > +*Pqq\qq > +ip q q, g \qq, £>+•••. (18) 

In target rest frame when the scattering energy W of the virtual photon is vary high 
and its virtuality Q 2 is rather large, we could know this from [TEj and [|nj that in the 
probability amplitude ip q q of decaying into a qq pair, there is a rather large component in 
which the longitudinal momentum fractions of the quark and antiquark are nearly equal and 
their transverse separation r t (r t (qq)) ~ ^(q)- After scattering off the target, its transverse 
expansion Ar± is very small, thus this qq pair could be treated approximatedly as a colour 
dipole in the interaction with nucleon (nucleus) target. In this case only those gluons which 
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comes from the target with transverse momentum k± ~ 0(Q) could probe and interact with 
this colour dipole, so now we can use perturbative QCD to discuss 7*-N(A) DIS process. 

Based on this physical observation, Levin et. al. || using the generalized Glauber 
(Mueller) pil formula for QCD diffractive coherent multi-scattering discussed how to get 
the gluon distribution function of nucleus, g A (x,Q 2 ) from that of nucleon, g N (x,Q 2 ). Since 
it is relevant for our later discussion, we will first give a schematically description about their 
approach. 

Since the PDF of nucleon are independence of the probe what we taken, in order to avoid 
complexity which comes from quark line loop in 7* — g fusion subprocess in 7*-N(A) DIS, 
Levin et al supposed a virtual "gluon probe" Q, which in target rest frame has very high 
energy and a rather large virtuality Q 2 . The probe Q decays into a gluon pair GG with 
nearly equal longitudinal momentum fraction and their transverse separation r± is rather 
small. The kinematics of this approach shown in Fig.[|. This pair interacts with a nucleon 
due to exchange of a gluon ladder diagram, where l± denotes transverse momentum of the 
gluon in the ladder attached to GG pair. The Bjorken variable x = where y/§ = W 
is the e.m. energy of incoming "gluon" Q. The change in r_|_ during the passage of GG 
pair through the nucleus Ar± oc R^jr, where E denotes the energy of the pair in target rest 
system and R is the size of target, the transverse momentum k± cx — . It is easy shown 
when W (and hence E) is very high and x <C one get Ar± <C r± and thus the condition 
for using Glauber (Mueller) approach in discussing the GG pair interaction with the nucleus 
is satisfied. Then as in || the gluon distribution function of nucleus is given by 

xG A (x, Q 2 ) = -,) d 4l d 4- / — 2{1 - e-^ G (*Vim) } . (19) 



7T 2 J X 1 J r\ J 7T 

x 



The term in curly brackets is the total cross section of the interaction of gluon pair GG with 
nucleus in the eikonal approach. In double logarithmic approximation (DLA) for perturba- 
tive QCD, the GG pair cross section with the nucleon ctqg(x, r ±) can be written (for N c = 3) 

m 



as 



3a, 



1 



4 



VGG(x,r ± ) = -^7T r ± (xg N (x, (20) 

4 TjJ 

The factor ^ in Eq.(^) comes from the wave function of the GG pair in probe Q and S(b±) 
is the nucleon profile function in nucleus (A) with impact parameter bj_ for Q-A eikonal 
scattering. Usually it takes Gaussian form: 
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Putting Eqs.(pO"D and ([Zl]) into Eq.(]T9|) and integrate over b± we obtain in DLA the gluon 
distribution of bound nucleon: 



x,Q 2 ) ee jxG A (x,Q 2 ) = ^J^rj ^{C + ln(k G (x',rl)) + E l (k G (x',rl))}, (22) 



xg 



ka{x, r]_) = ^ p2 ± xg N (x, — ) (23) 



where C is the Euler constant and E\ is the exponential integral and 

2>a s Aiir 2 

2R A , ± 

If Eq.(|22"D is expanded in small k G , the first term will correspond to the usual DGLAP 
equations in small x region (Born approximation of Mueller formula) while the other terms 
coming from GG pair coherently multi-scattering with nucleus what will take into account 
the shadowing corrections. 

For concrete comparison the xg A (x, Q 2 ) calculated, from Glauber (Mueller) approach 
with those from DQ 2 RM , we put Ql = 0.5 GeV 2 in Eq. (|22|) and calculate the xg A (x, Q 2 ) in 
2 < Q 2 < A GeV 2 , 10 -4 < x < 10~ 3 Q In this kinematic region we introduce the quantity: 

A ^ Max { -j^h^msMh^l h m 

10" 4 < a;<10~ 3 
2 < Q 2 < 4 GeV 2 

where xg A (x, Q 2 ) is calculated from Eq . (]22|) , and xg^Q 2RM (x,Q 2 ) is deduced from Eq.(|l|) 
and Table |, both of their input gluon distribution functions are LO of GRV parameterization 
sets |TT| . In addition, when we put the g A (x,Q 2 ) resulting from Eq. (p2j) into Eq.([I|), we get 



another Q 2 -rescaling parameter, £ A , for gluon distribution function from Glauber (Mueller) 
approach and we can compare it with those from Table |[ 



2 From the DLA condition which required for validity of Eq.(^), our Q 2 -range seems too low, 
but from the available data at such small x-region we found its Q 2 -range just so low or even 
lower. On the other hand from the fact of "precocious Q 2 -scaling" , which revealed experimentally 
what the Bjorken Q 2 -scaling in 7*-N DIS process arrived practically at much lower Q 2 values than 
those predicted by perturbative QCD, one may expected, as a first approximation, one could using 
Eq.p^) at such low Q 2 -range to calculate xg A (x,Q 2 ) 
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In Table HH we list the Af, and £* for some nuclei, from it we see the A; 4 values for 
■ — ' y 1 ^y 1 y 



all but H 1 and Be 9 nuclei^ are less than 5% and all the gluon Q 2 -rescaling parameters of 
DQ 2 RM , £g are all very close to those calculated from Glauber (Mueller) approach. 
In addition both the £r or £ A seems approach its saturation value, > 0.7 as A > 200, this 
is related to the fact that in 10~ 4 < x < 10~ 3 region the coherent length of the GG pair, 
Iqq ~ — « 10 3 — 10 4 fm are much large than radius of nuclei, Ra < 10 fm and Ra approch 
also its saturation value as A > 200, so the rescattering number of GG pair in nucleus also 
approach saturation. 

Here we want to point out that we could not expect in a similar simple manner as above 
formalism — Glauber (Mueller) approach in DLA — getting good results from calculating 
the sea quark distribution function of bound nucleon, xq$(x, Q 2 ). In other words, we could 
not take enough account the effect of nuclear shadowing on sea quark distribution functions 
from equations like Eq.(|T9|) since there are two points in this case: 

1. in small x- region, e.g. 10~ 4 < x < 10~ 3 sea quark distribution functions in nucleon are 
only | to | of that of the gluon. 

2. The relevant splitting function are P gq (x) = P 9 q(x) = | ) ; P gq (x) = 6[ ^_^ + 

+ x(l — x) + (y| — -jj)5(l — x)]. at small x limit their ratio approach |. 

These two factors would cause the corresponding power for sea quark in the exponential term 
in curly bracket smaller about one order of magnitude. Thus with this simple no mixing 
(between sea quark and gluon) approach the resulting effect of nuclear shadowing correction 
on sea quark is very weak and the resulting xqg(x, Q 2 ) would be incorrect. 

V. CONCLUSION AND COMMENTS 

We have adopted different approaches at different x-regions to demonstrate the 
DQ 2 RM for PDF of bound nucleon. For high x-region, 0.1 < x < 0.7 we used the 
distorted physical vacuum picture and the soliton bag model to justify the radius of a bound 



3 Notice Eq.(p2|) is not an evaluation equation but a formula for getting xG A (x,Q 2 ) of nucleus. 
In calculating xg A (x,Q 2 ) with Eq.(|22|), we have let the GRV parameterization set g N (x,Q 2 ) of 
nucleon as input. Since this input comes from phenomenology, it should have included effects of 
rescattering, thus in some sense, this calculating method have some double counting, but it is only 
important for nucleon and few light nuclei. 



11 



nucleoli always larger than that of free one. This effect causes the soften of PDF in bound 
nucleon, which turns to the <3 2 -rescaling behavior revealed phenomenologically on structure 
functions of bounding nucleon. We thus have demonstrated the validity of DQ 2 RM in high 
x-region. For low x-region < x < 10~ 3 according the formalism of Levin et al in 

Glauber (Mueller) coherently multi-scattering approach, we have shown the resultant gluon 
distribution function of bound nucleon, xg A (x, Q 2 ) agree very well with those from DQ 2 RM . 
Thus the validity for this phenomenological model are demonstrated. 

Of course, both the DQ 2 RM on its own and all our demonstrations about it are far from 
perfect. We want to make following observations about it. 

First, although our demonstrations for DQ 2 RM are physically reasonable, but it is far 
from satisfactory. We could not over the whole x-region using the perturbative QCD strictly 
and systematically justified this model, instead we adopted different approaches for different 
x-region. In particular we could not give a clear and physical argument for middle x-region, 
1CT 1 < x < 1CT 3 , where the contributions for valance quark, sea quark and gluon are 
comparatively in l-A DIS processes. 

Second, although the DQ 2 RM reflected approximatedly the objective reality of PDF in 
bound nucleon, but it is still a sketchy model. Indeed, if we have a detailed investigation on 
this model with perturbative QCD, it is possible to find that these Q 2 -rescaling parameters, 
£a( a = ViSi9)i are actually not constants, but perhaps like logarithmic, weakly depend 
both on x and Q 2 . 

Last, a very seriously problem is how to evaluate the non-perturbative QCD effects 
exerted on bound nucleons causing by the environment of nucleus. We have explained in 
Sec. Ill the EMC effect with the vacuum distored model in nuclei is just an example, but 
such non-perturbative QCD effects would be active over the whole x-region on PDF of bound 
nucleon. 

Obviously, to complete settle these problems from genetic QCD are very difficult. 
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TABLES 



TABLE I. The fitted values of rescaling parameters. 





C 12 


Ca 40 


Fe 56 


Sn 119 




1.03 


1.35 


1.14 


1.49 


is 


0.70 


0.67 


0.62 


0.60 


& 


0.86 


0.81 


0.76 


0.74 



TABLE II. £a and ta/tn for various nucleus. 



Nucleus 


u 


va/tn 


He 4 


1.223 


1.065 


Be 9 


1.242 


1.070 


C 12 


1.248 


1.072 


A 27 


1.299 


1.085 


Ca 40 


1.319 


1.091 


Fe 56 


1.352 


1.099 


Ag 107 


1.410 


1.114 



Au 197 1.696 " 1.180 



TABLE III. A list of A A , and if for various nucleus. 



A 


H 1 


Be 9 


C 12 


Al 27 


Ca 40 


Fe 56 


Cu 6 4 


Sn 119 


Ze 131 


Au 17 


p b 208 




12 


9 


4 


5 


4 


1 


2 


1 


1 


1 


1 




1 


0.9 


0.87 


0.84 


0.82 


0.78 


0.76 


0.74 


0.73 


0.72 


0.72 


^9 


1 




0.86 




0.81 


0.76 




0.74 
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FIG. 1. The ratio R A ' D (x, Q 2 ) for nuclei G 12 , Ga 40 , Fe 56 , Sn 119 versus the momentum fraction 
x with the corresponding experimental data [§,12]. The solid lines are the results of our model 
dashed lines are the results of the old Q 2 — rescaling model in ref. IT 
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FIG. 2. The ratio T A / D (x t ,Q 2 ) for nuclei C 12 , Ca 40 and Fe 56 versus x within region of 
0.025 < x t < 0.30, 4 < M /; - < 9GeV and = 40Ge^, where corresponding experimental data 

[13] are shown. The meaning of lines is the same as that in Fig.l. 



16 




17 




FIG. 4. The potential energy of c-field and the change of the bag constant due to the change 
of <7-field in the vacuum. 
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FIG. 5. Kinematics of the Glauber approach ||. 
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